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Abstract 

The asymptotic regimes of the TV-site complex Toda chain (CTC) 
with fixed ends related to the classical series of simple Lie algebras are 
classified. It is shown that the CTC models have much richer variety of 
asymptotic regimes than the real Toda chain (RTC). Besides asymp- 
totically free propagation (the only possible regime for the RTC), CTC 
allow bound state regimes, various intermediate regimes when one (or 
several) group(s) of particles form bound state(s), singular and de- 
generate solutions. These results can be used e.g., in describing the 
soliton interactions of the nonlinear Schrodinger equation. Explicit 
expressions for the solutions in terms of minimal sets of scattering 
data are proposed for all classical series B r - D r . 
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1 Introduction 



The Toda chain model |, |, |, |, |, |, [7|, |, |, 

d 2 q k 



dt 2 



exp (g fe+1 - q k ) - exp (g fc - q k -i), (1) 



is one of the paradigms of integrable nonlinear chains and lattices. It has 
been thoroughly studied for a number of initial and boundary conditions, 
such as: 

• fixed ends boundary conditions, i.e. go = — Qn+i — °°; this will be the 
case we are interested in; 

• infinite chain — oo < k < oo with lim^-oo Qh = and lim^oo^ = 
const, or equivalently, lim fc ^- too (g fc+1 — q k ) = 0; 

• quasi-periodic boundary conditions q k +N = Qk + c, where c = const. 

This model appeared first in describing the oscillations of a one-dimensional 
crystalline lattice J/J. Since then many other applications became known, see 
e.g. !• 

The model ([!]) is directly related to the algebra sl(N), where N is the 
number of sites of the chain. Most of the references cited above are devoted to 
the case when q k (t) are real- valued functions. That is why for definiteness we 
will call this model real Toda chain (RTC). Other generalizations of the RTC 
are related to: a) simple Lie algebras; b) affine (or Kac-Moody) algebras; 
c) 2- dimensional generalizations. 

Another possibility for generalizing the RTC, which as far as we know has 
not been investigated, is the complex Toda chain (CTC). We see two main 
reasons for this: 

(i) the generic solutions of the CTC are readily obtained from the ones of 
RTC by simply making all parameters complex. In fact, technically solving 
the RTC requires additional efforts in ensuring that the scattering data of 
the Lax matrix L are real- valued. 

(ii) the CTC has not been known to have physical applications. 



Recently however, it was discovered [O, IT4L [L5| that the CTC describes 



the iV-soliton train interaction of the nonlinear Schrodinger (NLS) equation 
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in the adiabatic approximation. The variables q k (t) in (|I]) are related to the 
solitons parameters by 

q k = -2v £ k + i (2/i £ fc - S k ) + const (2) 

where and 5 k characterize the center-of-mass position and the phase of the 
k-th soliton in the chain; vq and /j,q are the average amplitude and velocity of 
the soliton train. Such soliton trains and their asymptotic behavior appear 
to be important for the needs of soliton based fiber optics communications. 

In what follows we will view the CTC as a model of "complex" particles 
characterized by a coordinate and phase which, in analogy with (g), are 
related to Req k (t) and lmq k (t) respectively. 

Another reason for the present paper is in the fact, that along with the 
similarities between the solutions of CTC and RTC, there are also important 
qualitative differences between the asymptotic properties of these solutions. 

The purpose of the present paper is to derive analytically the large- time 
asymptotics of the solutions to the CTC models: 

§^t^=T.H,e-'^, (3) 

related to the classical series of simple Lie algebras extending the results of 
[fL6| , [1T|| . Here n g = {cti, . . . , a r } is the set of simple roots of the algebra g, 
q(t) is a complex-valued function of t taking values in the Cartan subalgebra 
h of g; H k form a basis in h dual to the orthonormal basis {e k }' k=1 in the 
root space E r and r is the rank of g; H a = X^ =1 (a, e k )H k . For more details 
about the structure of the simple Lie algebras, see e.g. |P| . Equation ([[]) is 
a particular case of (f|) for the A r series, i.e., for g ~ sl(r + 1). These results, 
and especially the ones for the series A r , can be used as a tool for deriving 
the asymptotic behavior of the iV-soliton trains from the initial set of soliton 



parameters |T3|, [L4|, |15| . 

We also specify the minimal sets of scattering data T g for L which deter- 
mine uniquely both L and the solutions of @ and obtain explicit expressions 
for the solutions of (|3|) in terms of T g . 



2 Comparison between RTC and CTC 

Since the paper of Moser || on the finite nonperiodic real Toda lattice, there 
have been proposed many methods for solving ([!]) for the various choices 
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of the initial and boundary conditions, see |2], [| |, ||, |6|, 0, ^, |, [H], |TT 
and the numerous references therein. The RTC model was also extended 



to indefinite metric spaces [I5| and was shown to possess singular solutions 
"blowing up" for finite values of t. These models can be viewed as special 
cases of the CTC in which part of the a k s defined below are real while the 
rest are purely imaginary. 

As we already mentioned, a number of properties of the CTC are obtained 
from the corresponding ones of RTC trivially by assuming the corresponding 
variables complex. We list the four most important of them below: 

51) The Lax representation: There are several equivalent formulations of 
the Lax representation for @. Below we will use the "symmetric" one: 

r 

L(t) = Y,(hH k + a k (E ak +E_ ak )), (4) 

fc=l 

M{t)= l -Y,a k {E ak -E_ ak ) 1 (5) 
A k=i 

where a k = \ exp (— (q, a k )/2) and b k = \dq k jdt\ for g ~ sl(N) we have 
a k = | exp ((qk+i — g&)/2). It is well known that to each root a G A g C E r 
one can put into correspondence the element H a e h. Analogously, to q(t) = 
Reg(t) + ilmq(t) there corresponds the vector q(t) = Req(t) + ilmq(t), 
whose real and imaginary parts are vectors inJ^ r . 

52) The integrals of motion in involution are provided by the eigenvalues, 
Cfc, of L. 

53) The solutions of both the CTC and the RTC are determined by the 
scattering data for L = L(0). When the spectrum of L is nondegenerate, 
i.e. Cfc 7^ Cj f° r 7^ 3, then this scattering data consists of 

T = {Ci,---,CN,ri,...,r N }. (6) 

where r k are the first components of the corresponding eigenvectors r k = V\ k 
of L in the typical representation R{lo\) of g, N = dimi?(o;i): 

L V = VZ, Z = diag(Ci,...,Cv). (7) 
They are determined (up to an overall sign) by the normalization conditions: 

E(Kfc) 2 = (^ (s) V s) ) = l, s = l,...,N, (8) 
fc=i 
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see |, |, |; then V T = V' 1 . 

S4) Lastly, the eigenvalues of L uniquely determine the asymptotic be- 
havior of the solutions; these eigenvalues can be calculated directly from the 
initial conditions. We will extensively use this fact for the description of the 
different types of asymptotic behavior. 

However, there are important differences between the RTC and CTC, 
especially the asymptotic behavior of their solutions. Indeed, for the RTC, 
one has [[3], || that, both the eigenvalues, ( k , and the coefficients, r k , are 
always real- valued. Moreover, one can prove that ( k ^ Cj f° r k ^ j, i.e. 
no two eigenvalues can be exactly the same. As a direct consequence of 
this, it follows that the only possible asymptotic behavior in the RTC is an 
asymptotically separating, free motion of the particles (solitons). 

This situation is different for the CTC. Now the eigenvalues ( k = K k + ir] k , 
as well as the coefficients r k become complex. Furthermore, the argument of 
Moser can not be applied so one can have multiple eigenvalues. The collection 
of eigenvalues, £fcj still determines the asymptotic behavior of the solutions. 
In particular, it is K k that determines the asymptotic velocity of the fc-th 
particle (soliton). For simplicity, we assume ( k 7^ Cj f° r k ^ j. However, this 
condition does not necessarily mean that K k ^ Kj. We also assume that the 
/tfc's are ordered as: 

^i < ^2 < • • • < Kn- (9) 

Once this is done, then for the corresponding set of N particles (train of N 
solitons), there are three possible general configurations: 

Dl) K k 7^ Kj for k 7^ j. Since the asymptotic velocities are all different, 
one has the asymptotically separating, free particles (solitons). 

D2) K\ — K2 — ■ ■ ■ — k>n- In this case, all N particles (solitons) will move 
with the same mean asymptotic velocity, and therefore will form a "bound 
state". The key question now will be the nature of the internal motions in 
such a bound state. 

D3) One may have also a variety of intermediate situations when only one 
group (or several groups) of particles move with the same mean asymptotic 
velocity; then they would form one (or several) bound state (s) and the rest 
of the particles will have free asymptotic motion. 

Obviously the cases D2) and D3) have no analogues in the RTC and 
physically are qualitatively different from Dl). The same is also true for the 
special degenerate cases, where two or more of the C^'s may become equal 
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and for the singular solutions. These cases will be considered briefly below. 



3 Solutions of the CTC 

The solutions for the CTC can be obtained formally from the well known 
ones for RTC by inserting the corresponding complex parameters. We remind 
first the solution of the RTC for g ~ sl(N), see |2D| and the references 
therein. Here we fix up the mass center at the origin by 

N 

= o- (io) 

k=l 

The velocity of the center of mass is given by tr L = Y, k =i Cfc = 0, due to 
L G sl(N). Then the solution has the form: 

?*(f)=fc(0)+ln^- (11) 

where A = 1, 

N 

Mt) = T,rle- 2 ^ (12) 
fc=i 

Mt) = E (r il r l2 ...r i j¥(I 1) I 2) ... ) ye- 2 (V-+y (13) 

l<«i</ 2 <-</fc<A r 

and 

N 

A N = W 2 (l, 2, . . . , N) J] r 2 k = exp (-iV gi (0)). (14) 

fc=i 

By W(h, . . . ,lk) we denote the Wandermonde determinant: 

W(h,...,l k )= J] (2C.-2Cp)- (15) 

s > p 

s,p e {h, . . . ,4} 

The solutions of (0) for real- valued q(t) are well known in several different 
formulations. We remind here the formula, which effectively is contained in 

I: 

(q\t),u k ) - (q(0),uj k ) = \n(u k \exp(-2L t)\u k ), (16) 
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where Uk are the fundamental weights of g. The fact that the large time 
asymptotics of quit) for the g ~ sl(N) RTC have the form: 

lirn {q k {t) - v±t) = P± (17) 

t — >±oo 

where the asymptotic velocities v£ = —2(k and v k = —2C,N+i-k have been 
derived by Moser ||. He also evaluated the differences (3 k — P^-k+i which 
characterizes the interaction in the RTC model. 

The minimal set of scattering data for g ~ sl(N) is obtained from (0) by 
imposing on T the restrictions J2k=i Cfc = and 

N 

£r| = l. (18) 

fe=i 

Note that exp (— qi(0)) is expressed through T by (0). 

For the RTC related to the other classical series of Lie algebras it is known 



TK 137711 that 



lirn (q(t)-v ± t)=(3 ± , (19) 

t^±oo 

v ± G h, /? ± G h and t> + = u>o(f~), where w is the element of the Weyl group, 
which maps the highest weight of each irreducible representation of g into 
the lowest weight. The action of Wq on the simple roots is well known ||18||: 



w (a k ) = -a~ k , (20) 

k = r — k + 1 for A r , k = k, k = 1, . . . , r for B r , C r and, when r is even, also 
for D r . When g ~ D r and r is odd k = k for < r — 2, and tuo(av-i) = — «r, 

W («r) = -CKr-1- 

What we will do below is to: i) specify how minimal sets of scattering data 
T 9 can be extracted from @; ii) find explicit expressions for /3 k = (/5 ± ,e/ c ) 
for each of the classical Lie algebras in terms of T g . 

Like in the sl(N)-c&se ( k and r k are the eigenvalues and the first compo- 
nents of the eigenvectors of L in the typical representation, namely: 

L V = VJ2 ( k H k , (21) 
fc=i 

The requirement that L (and as consequence, L) belongs to one of the alge- 
bras in the B r or C r series imposes on q k the following natural restrictions: 

Qk = —qN-k+i (22) 
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which leads to 

0-k = dN-k (23) 

b k = -b N +i-k- (24) 

Thus the solutions for g ~ B f . and C r can formally be obtained from the 
ones for sl(N) (|TT|)-([H]) by imposing on them the involutions (^)-(|23|). So 
we have to find out what are the restrictions on T imposed by (p2j)-(|24|); 
this will provide us with the corresponding minimal set of scattering data 
T g , which must obviously contain only 2r parameters. It is easy to find that 
= —(k, so only r of them are independent. It is not so trivial to derive the 
corresponding relations which reduce the number of the coefficients r^. Our 
analysis shows that: 

r k n = exp (-<7i(0))wfc, k = l,...,r (25) 

where k = N + 1 — k and Wk is expressed in terms of Below we provide 
the explicit formulas for Wk for each of the classical series B r , C r and D r . 

B r -series: N — 2r + 1. Note that in this case ( r+ i = 0, and in addition 
to©: 

1 



r 



2 



> + i-exp(- gi (0)) 22r(CiC2 (26) 
and the expression for Wk is provided by: 

^ fe— i 1 r 1 

Wk = W k S <F^a J? +1 icF^S' (27) 

Inserting (]25])-(|2~7|) into (0) we obtain a quadratic equation for exp (— <7i(0)), 
so it can be expressed in terms of T g . 
C r -series: N = 2r. Here 

w ' = -^H4^,n +1 4cFic?' (28) 

Like for the B r -series exp (—51(0)) is determined from fll~8l). 
D r -series: iV = 2r. Here 

m * = Hi^ n +1 icFicf- (29 » 
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Again exp (— qi(0)) is determined from (|T8| ) and (p5|). The derivation of the 
solution for this series requires some additional efforts. The main problem 
here is to find explicit parametrization for the right hand sides of ( |IED for 
k = r — 1 and r in terms of r^, which characterize the matrix V in ([H]) in 
the typical representation. Skipping the details we just present the result, 



namely that q^ with k = 1, . . . , r — 1 are again given by (|lT|)-(|r3|) where 
Cfc = — Ck an d r k are restricted by (^), (|^). For k = r the solution for q r (t) 
is provided by flTT|) with 



l<h<l 2 <...<lr<N 

and 



fh = 5(1+ ■ (3D 

The proof of these formulas is based on detailed analysis of the properties of 
the fundamental representations R(u>k) of g and of their tensor products. 

We remind that these formulas are valid both for RTC and CTC In the 
latter case one should be careful to avoid the subset of singular cases, when 
one or more of the functions A^it) may develop zeroes for finite values of t, 
see also the discussion in Section V below. 



4 Large time asymptotics 

Let us now express large time asymptotics of qk{t) in terms of the minimal set 
of scattering data T g and analyze the different types of asymptotic regimes. 

As we mentioned above, we view the CTC as a model, describing non- 
trivial scattering of iV "complex particles" so that Regfc(t) and \mqk{t) cor- 
respond to their coordinates and "phases" . 

Dl) Asymptotically free states 7^ k,j for k 7^ j. It is specified by 
imposing on (k = + iVk, k — 1, . . . , N the so-called sorting condition: 

Ki < K 2 < . . . < kn- (32) 

Now we have to express /3 k in terms of T g . Skipping the details we list 
the results for the classical series of simple Lie algebras A r - D r . 
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A r -series: iV — r + 1. The corresponding T s is formed from {^,^^^=1 
by imposing Y,k=i Cfc = and the normalization condition ([Tj|). Then 

A;— 1 

/3+ = 9l (0) + lnr^ + In J] (2& - 2G) 2 (33) 

JV 

= 9l (0) + In rf + In (2£ - 2Cs) 2 , (34) 

s=fc+l 

where /c = N + 1 — k. 

Now it is easy to calculate the shift of the relative position which is the 
effect of the particle interaction. Naturally these shifts are also complex- 
valued. If qk correspond to NLS solitons then (3£ — (3 k will describe the 
shifts of both the relative positions and phases of the solitons, see formula 
(H). Note that in the class of regular solutions the particles do not collide, 
i.e. their trajectories do not intersect. Since we have ordered the particles by 
their velocities assuming K\ < k 2 < . . . < kn, so for t — > 00 the k-th particle 
will move with velocity —2ku- For t — >• — 00 however we find that due to 
the interaction now the k-th. particle moves with velocity — 2KN-k+i = — 2/^. 
This is the so-called "sorting property" characteristic for the RTC. If we 
identify the k-th particle (soliton) by its velocity then its shift of position 
will be given by the real part of the relation: 

/^-/# = £^ ln (2C;-2G) 2 , (35) 

where €jk = 1 for j > k and €jk = — 1 for j < k. The imaginary part of ( jffi ) 
will provide the shift in the phase of the corresponding particle (soliton). 
Equation ( |35|) is a natural generalization of the corresponding result of Moser 
[|] for RTC. 

The asymptotics for the B r , C r and D r series have the form ([19D : 

lim (q k (t)±2C k t)=/3t. (36) 

t — >±oo 

with 

fc— 1 

ft = <Zi(0) + \nr 2 k + In I](2C - 2( k ) 2 (37) 
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w 2 k ~ l 

Pk = -9i(0) + In -i + In J] (2C - 2C*) 2 (38) 

r k s=l 



for /c = 1, . . . , r. 



The only exception is for the series D r in the case of odd r, k = r and 
t — > — oo which reads: 



lim (q r (t) + 2( r t) = (3; (39) 

t — > — OO 

where ^ 

/3- = 9l (0) + lnr r 2 + In n(2Cr - 2Cs) 2 . (40) 

s=l 

D2) Bound states: /cj = k,2 = ■ ■ ■ = k>n — 0. Now all particles 
will move with the same mean asymptotic velocity for both t — ► oo and 
t — > — oo; by Galilean transformation this velocity can always be made 0. 
The individual velocities of the particles oscillate around the common mean 
value. In other words we find that all N particles generically do not separate 
but form a bound state with a large number of degrees of freedom. The 
explicit solutions for now do not simplify even for t — > ±oo. Nevertheless 
two features are worth noting. 

The solutions will be periodic functions of t provided the ratios (% — 
Vm)/{Vk — Vj) are rational numbers for all k, m and j. In some cases they 
can also become singular, see the next section. 

Important for possible physical applications is the so-called quasi-equidistant 
regime in which the distances between the neighboring particles, i.e. Re (qk+i{t) — 
quit)) oscillate with very small amplitude; with rather good accuracy they 
could be considered constant. The fact that such regimes are possible and 
describe adequately the behavior of certain NLS soliton trains is shown in 

0- 

D3) Mixed regimes. As we mentioned above there are a number of 
intermediate cases. Here we start with the case when m + 1 out of the 
particles form a bound state, i.e.: 

Kl < . . . < K k = . . . = K k+m < Kk+m+1 < ■ ■ ■ < K N (41) 

and rji ^ rjj for j G {k, k + 1, ... k + m}. Skipping the details we present 
the results for the case g ~ sl(N) and m = 1, i.e., only two of the particles 
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form a bound state. Particles with numbers different from k, k + 1 are free, 
and for k-th and k + 1-st have 

q k+a (t) = 9l (0) + u fc t + /3 fe + +a (0 + 0(e-"), (42) 

for i — > oo and 

g^_ fc+a (t) = 9l (0) + u fc t + + O(e^), (43) 

for t -> -oo. Here a = 0, 1, w fc = + Cfe+i)> 

z/ = min(/t fc - /c fe _i, K fe+2 - Kfc+i), 

and 

Pk + a = Pk + (0 + aln(2C fc - 2C fc+ i) 2 , (44) 

B±(t) = In (2 cos(( % - Vk+1 )t - i 7 ±)) , (45) 

/3± = In I r fc r fc+1 JJ^C, - 2G)(2& - 2 Cfe+1 ) (46) 
where EL/ = ifo 1 and n.,*" = Utk+2- 

5 On the singular and degenerate solutions 
of CTC 

It is only natural that some of the CTC solutions do not enjoy all the good 
properties of the RTC ones. We mentioned already two of them: 

PI) The elements of T g for the RTC are real- valued; 

P2) The eigenvalues ( k are pairwise different, i.e. ( k ^ Q for k ^ j. 

An immediate consequence of PI) and P2) is the fact that the solutions 
q k (t) of the RTC are regular functions for all finite values of t. 

Generalizing to CTC we loose both properties PI) and P2). As a result, 
besides the regular solutions, we obtain also singular and degenerate solu- 
tions. Obviously, all solutions leading to the Dl regime are regular. Even if 
we assume that property P2) holds, we may still have singular solutions. 
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Indeed, from equations ( fj5|) - (H7[) we see that in the "oscillating part" 
of the motion Bf(t) is periodic. If in addition the parameters in T g are such 
that Re 7^ = in (^) then B k {t) will develop singularities for finite values 
of t. The same holds true for the functions qk{t)'- there exist submanifolds of 
T g for which qu{t) become singular for finite values of t. This fact has been 



noted in |19| for RTC on spaces with indefinite metric and (or) with purely 
imaginary interaction constant. Such RTC can be viewed as particular case 
of the CTC when one or several of the functions ak{t) are purely imaginary, 
while the others remain real. 

Let us now examine the degeneration, i. e. the case when the property 
P2) is violated and some of the eigenvalues of L become equal. In || |22| it 



is proved that the spectrum of Lq is simple for real s/(iV)-Toda chain. We 
state the following generalization: 

Lemma 1 Let us consider Lq for the complex Toda models related to the 
classical series of simple Lie algebras A r , B r , C r , D r . If Lq does not contain 
Jordan cells, then its spectrum is simple. 

Therefore the degeneration can take place only for CTC models and only 
provided in the diagonalization of L Jordan cells are present. Let us have 
g ~ sl(N) and let Lq has a 2 x 2 Jordan cell, Ci = C2 = C Then Lq has an 
eigenvector t>W and an adjoint eigenvector v^: 

L oV W = (48) 

L v^=Cv^+v^\ (49) 

where can be expressed as linear combination of v^(() and its first 
derivative with respect to (. The corresponding A k (t) besides the standard 
exponential terms will contain also terms of the form te~ 2C - t . More generally, 
if the degeneracy is of higher order, i.e., Ci = • • • = Cm then we will need 
linear combinations of v ^(C) and its derivatives with respect to ( of order 1, 
. . . , m— 1 and will contain terms of the form t p e~ 2n ^ with p = 1, . . . , m — 1 
and n = 1, . . . , min(fc, m), see also 



In particular, if we have complete degeneracy (i.e., all (k are equal and 
equal to zero) the solution of the s/(iV)-CTC is expressed through Ak, which 
are polynomials of degree k(N — k). Their coefficients depend 011 iV - 1 
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constants fk, k = 1, . . . , N — 1. For example, for N = 3 and Ci = C2 = C3 = 
we get: 

A 1 (t) = -^t 2 + f 1 t + f 2 , (50) 

A 2 (t) = -\t 2 + ht - fl - f 2 . (51) 

and A3 = 1. Obviously, these solutions will be regular if have complex 
roots and will develop singularities if one (or more) of their roots are real. 



6 Conclusions 

We have shown that the iV-site CTC with fixed ends has richer variety of 
asymptotical regimes than the RTC. 

In particular, we showed that the CTC allows solutions in which the 
"complex" particles form regular bound states. This could be important 



for the applications to soliton interactions, see e.g. 0] where this is done 
for N = 3. Then the problem to determine the initial soliton parameters 
(i.e., the values of and bk) is reduced to the analysis of the characteristic 
equation for L : 

N 

det(L -C) = E^C fc (52) 

k=0 

Pn = 1, Pn-i — and to the requirement that the eigenvalues of L be 
purely imaginary. This is an algebraic problem which often can be solved 
analytically. It allows one to determine the set of initial soliton parameters 
in such a way, that the solitons will not only form a stable bound state, 
but also will propagate quasi- equidistantly. This type of propagation is of 
importance for soliton based fiber optics communications f21 ]. 



We have studied also the asymptotic regimes for the CTC related to the 
other simple Lie algebras. We proposed explicit solutions to these models in 
terms of the minimal sets of scattering data T g . The degenerate and singular 
solutions of the CTC which have no counterparts in RTC are also briefly 
analyzed. 
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